#### Figure S9: Bivariate Support for Courts ####

# Libraries
# library(here)
# library(rio)
# library(tidyverse)
# library(ggthemes)
# library(lubridate)
# library(estimatr)

# data_pnas = import(here("Data","data_pnas.rds"))
# source(here("Code","Functions","funcs_theme.R"))
# source(here("Code","Functions","funcs_bvplot.R"))

fig_s9 = plot_bv_maga(data_pnas, norm_judgesc)

print(fig_s9)

ggsave(here("Plots","Supplementary","figure_s9.png"),
       fig_s9,
       units = "in",
       width = 6, height = 6,
       dpi = 600)